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ABSTRACT 

Recent ROSAT surveys suggest that galaxies can constitute the new class of faint 
sources required to explain the full phenomenology of the cosmic X-Ray Background 
(XRB). To test this hypothesis without resorting to optical identifications, we compute 
the two-point cross-correlation function (CCF) estimator W xg (9) between 62 Einstein- 
IPC fields (.81-3.5 keV) and the APM Northern galaxy catalog (13.5 < E < 19.0). 
At zero-lag (8 = 0), we detect a 3.5<r correlation signal with an amplitude of 
W xg (0) = .045 ± .013. This signal passes a series of control tests. At non-zero lag 
(6 > 0), the angular dependence of W xg has two main features: the main signal for 
9 ^ 4', and an almost flat plateau with an amplitude of W xg (6 ^ 4') ~ .015. When 
fields with galaxy clusters as Einstein targets are removed, the plateau virtually 
disappears, and the zero-lag amplitude becomes W xg (0) = .029 ± .013. We develop a 
simple, 2-dimensional formalism to interpret the CCF which takes into account the 
point-spread function of the imaging X-ray detector. Three distinct effects can produce 
a correlation signal: the X-ray emission from galaxy themselves, the clustering of 
galaxies with discrete X-ray sources, and the clustering of galaxies with diffuse X-ray 
emission. It is likely that the plateau at large angles is due to the last effect through 
the residual diffuse X-ray emission from clusters of galaxies. We do not detect any 
significant clustering between discrete X-ray sources and galaxies. Using only the fields 
with non-cluster targets, we find that the mean X-ray intensity of APM galaxies in the 
.81-3.5 keV band is (2.2 ± 1.1) x 10~ 6 cts s" 1 arcmin , corresponding to 1.5 ± .8% 
of the XRB intensity. The mean X-ray flux of galaxies with (E) = 17.5 ± .3 is then 
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(8.1 ± 4.7) x 10~ 16 ergs s -1 cm -2 . This agrees within la with the X-ray flux expected 
from earlier direct studies of brighter, nearby galaxies, which were shown to result in a 
total integrated galaxy contribution to the XRB of about 13%. We discuss how this 
powerful cross-correlation method can be used to measure the flux of X-ray sources 
well below the detection limit of X-ray instruments, and, perhaps, to probe otherwise 
undetectable faint diffuse X-ray emission. 

Subject headings: Cosmology: Diffuse Radiation - X-Rays: Galaxies - Galaxies: 
Statistics - X-Rays: General - Methods: Statistical 

1. Introduction 

In spite of significant recent progress, the nature of the cosmic X-Ray Background (XRB) 
remains one of the outstanding puzzles in astrophysics (for reviews, see Fabian & Barcons 1992; 
Zamorani 1995; De Zotti et al. 1995; Hasinger 1996). The observation of the distortion of 
the Cosmic Microwave Background spectrum with COBE (Mather et al. 1990) precludes the 
possibility that the hard XRB (e ^ 2 keV) originates from a homogeneous diffuse hot plasma. The 
alternative explanation, namely a superposition of discrete sources, is supported, at least in the 
soft band (e <J 2 keV), by deep ROSAT surveys (Hasinger et al. 1993; Vikhlinin et al. 1995a,b,c) 
in which about 60% of the XRB is resolved into discrete sources, of which, in turn, about 60% 
are AGN. However, reconciling the full phenomenology of the soft and hard XRB with known 
properties of AGN has led to several difficulties (Hasinger 1996). First, no single class of known 
objects has a spectrum consistent with that of the XRB. Also, the surface density of faint ROSAT 
sources is higher (Hasinger et al. 1993) than expected from the luminosity function of known AGN 
(Boyle et al. 1993). In addition, fluctuation analyses with GINGA and HEAOl (Butcher et al. 
1994; Mushotzky & Jahoda 1992) imply more sources in the hard band than expected from the 
ROSAT deep counts by a factor of 2-3. 

New detailed unified AGN models stimulated by recent observations could provide a 
self-consistent solution (for a review, see Setti & Comastri 1996). Another possibility is that a 
new class of sources becomes important at low fluxes. It has been recently proposed that galaxies, 
whose fraction is observed to rise at faint fluxes in optical identifications of deep ROSAT surveys, 
could constitute this new class (Jones et al. 1995; Carballo et al. 1995; Boyle et al. 1995; 
Romero-Colmenero et al. 1996; see also Hasinger 1996). 

One way to study the contribution of galaxies to the XRB in addition to the time-consuming 
optical-identifications, is to cross-correlate the XRB intensity with galaxy catalogs. A sensitive 
statistical measure of the correlation at a given angular lag 9 is provided by the two-point 
cross-correlation function (CCF) W xg (6). The 9 = value of W xg , or "zero-lag" value, is often 
used, since, in most cases, most of the signal resides at small angles. However, the full correlation 
signal is provided by the "non-zero lag" (9 > 0) sector of W xg . 
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This approach was initiated by Jahoda et al. (1991) who computed W xg at zero lag between 
2-10 keV HEAO 1 A- 2 maps of the XRB and the UGC and ESO galaxy catalogs. They found a 
significant correlation signal with an amplitude of W xg (0) = (3 ± 1) x 10~ 3 , and concluded that as 
much as 50-70% of the XRB could be due to X-ray sources associated with present-epoch galaxies. 
Lahav et al. (1993) confirmed the correlation signal by applying the technique to 4-12keV GINGA 
maps of the XRB compared with the UGC and IRAS galaxy catalogs. However, they pointed out 
that a correct interpretation of the correlation requires taking source clustering into account. Their 
revised estimate of the fraction of the XRB due to a non-evolving population of X-ray sources 
spatially associated with local galaxies is 30 ± 15%. From a careful modelling of the correlation 
found between the HEAO 1 A-2 maps and the IRAS galaxy catalog, Miyaji et al. (1994) derived a 
local X-ray emissivity which correspond to about 20% of the XRB without evolution and which 
can be accounted for by AGN alone (see also Barcons et al. 1995). A similar fraction (10-30%) 
was obtained by Carrera et al. ( 1995a, b) from a measurement W xg {6) at non-zero lag, between 
GINGA and several galaxy catalogs. 

The aforementioned correlation studies were focused on the hard (E > 2 keV) XRB. Because 
only data from non-imaging X-ray instruments were available in this band, this work was restricted 
to rather large angular scales {6 > 1 deg). In addition, only catalogs of bright, nearby galaxies were 
employed. More recently, Roche et al. (1995) detected a correlation signal between ROSAT-PSPC 
fields and faint galaxies (B < 23). The high angular resolution of ROSAT allowed them to probe 
small angular scales (~ 15"). However, their consideration of only three ROSAT fields limited their 
statistics. By applying a detailed correlation formalism to the ROSAT results, Treyer &; Lahav 
(1995) concluded that about 22% of the .5-2 keV XRB is due to faint galaxies with 18 < B < 23. 

In order to test these results at intermediate angular scales (~ 1') and with improved 
statistics, we have cross-correlated 62 Einstein-IPC fields with the APM Northern Galaxy catalog 
(13.5 < E < 19). To extract the full correlation information, we have computed both the zero and 
non-zero lag CCF. The improved statistics allows us to give a better estimate of the statistical 
uncertainty in the correlation signal. Several control experiments provide added confidence in the 
significance of our results. The Einstein energy band (.1-3.5 keV) is sufficiently broad for testing 
the energy dependence of the signal. In addition, we develop a simple, 2-dimensional formalism 
to interpret the angular dependence of W xg which takes into account the Point Spread Function 
(PSF) of the imaging X-ray instrument. 

This paper is organized as follows. We first describe the X-ray and galaxy data (§§). We then 
outline our procedure for measuring the correlation function estimator W xg (§||). In §||, we present 
the zero and non-zero lag results, along with the analysis of the various control samples. We then 
present our interpretation formalism (§||) and discuss the implications of the detected correlation 
for the contribution of galaxies to the XRB (§^). The details of the calculations involved in 
the interpretation of W xg are relegated to appendix |A[ In appendix [B|, we explicitly evaluate 
two multi-dimensional integrals which appear in this interpretation. Finally, a treatment of the 
statistics of W xg is presented in appendix The zero-lag results have already been reported in a 
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brief form in Refregier et al. (1995). 

2. Data 
2.1. X-Ray 

The original X-ray data consisted of northern (Dec> —5° ) Einstein Imaging Proportional 
Counter (IPC) fields (Giaconni et al. 1979) at high galactic latitude (|6| > 25°). We originally 
selected 77 random fields with typical exposure times greater than 5000 s. We avoided fields 
containing bright sources and/or large areas of diffuse emission. The mean exposure time was 
(t exp ) ~ 1.4 x 10 4 s. In order to focus on the unresolved component of the XRB, sources were excised 
down to 3.5<7. This corresponds to a mean excision flux of (S exc (.81 — 3.5keV)) = (4.8 ± .5) x 10~ 14 
ergs cm -2 s _1 . The source excision and flat-fielding procedures are described in Hamilton et al. 
(1991) and Wu et al. (1991). The type of target was recorded for each field and classified into 
four broad categories (see Table ||): galactic sources, extragalactic discrete sources (galaxies, QSO, 
etc), clusters (and superclusters) of galaxies, and deep fields. This allowed us, in particular, to 
study subsamples of fields with and without clusters of galaxies as targets. 

The photons in each one-square-degree field where binned into 64" pixels. The Point-Spread 
Function (PSF) for the IPC is close to a gaussian with standard deviation a ps f ~ .5'. A typical 
X-ray field is shown as the grey-scale map in figure |l[ The blank region in the 4'-wide criss-cross 
pattern, due to the window support ribs, has been excluded from our analysis. We considered two 
energy bands: "soft" (.16-. 81 keV) and "hard" (.81-3.50 keV). The hard band is to be considered 
as the data set of interest, whereas the soft band, which is dominated by emission from the Milky 
Way is used as a control data set. In general, fields in the Einstein archive have different levels 
of contamination from solar X-rays scattered into the telescope by the residual atmosphere. This 
is parametrized by the Viewing Geometry (VG) flag with VG=1 (VG=3), corresponding to low 
(high) contamination. For the real data set, only exposures with relatively low solar contamination 
(VG=l-2) were considered. Another control data set was constructed by keeping higher solar 
contaminations (VG=2-3). Even though the low and high contamination data sets both contain 
VG=2 exposures and thus overlap, they can be compared to test the qualitative dependence of 
the correlation signal on solar contamination. 

For the Einstein-IPC, counts result from several effects (Wu et al. 1991). The instrumental 
background due to cosmic-ray particle contamination and calibration source leakage amounts to 
Impart _^_^cahb^ _ ^,1 -j- .1) x 10 -4 cts s _1 arcmin -2 . The latter and following intensities correspond 
to the hard band (.81-3.5 keV) and to the subset of our fields with targets other than galaxy 
clusters. The solar contamination for VG = 2 — 3 exposures is well described by a bremmstrahlung 
spectrum with kT = .25 keV and produces, on average, about .19 x 10~ 4 cts s -1 arcmin -2 (Wang 
et al. 1991). We thus estimate that, for VG = 1 — 2 exposures, this contribution produces 
(i so1 ) = (.1 ± .1) x 10 -4 cts s" 1 arcmin . The XRB, including the contribution from the Milky 
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Way, produces (i XRB ) = (1.4 ± .1) x 10 -4 cts s -1 arcmin -2 . By integrating the number-flux 
relation measured by Gioia et al. (1990) in the context of the Einstein Medium Sensitivity Survey 
(EMSS), we estimate that the source excision mentioned above amounts to the removal of an 
average of (i exc ) = (.10 ± .03) x 10 -4 cts s -1 arcmin -2 . (We have assumed a mean power law 
spectrum with index .7 for the sources and Nh = 3 x 10 20 cm -2 to convert from the EMSS band). 
The final expected count rate is thus (i) = {i part + i calih + i XRB - i exc ) = (2.5 ± .3) x 10 -4 cts 
s -1 arcmin -2 . For our data set, we measure a mean intensity of (i) = (2.73 ± .01) x 10 -4 cts s -1 
arcmin -2 , consistent with the expected count rate. 

2.2. Galaxies 

The APM Northern Catalog (Irwin & McMahon 1992; Irwin et al. 1994) consists of scans 
of E and O Palomar Observatory Sky Survey (POSS-I) plates by the SERC Automatic Plate 
Measuring (APM) machine in Cambridge (Kibblewhite et al. 1984). Objects are assigned position 
coordinates, E and O apparent magnitudes, a morphological class and various morphological 
parameters. The objects in the raw APM object catalog were classified on the basis of a 
compactness parameter similar to that used by Maddox et al. (1990b). The morphological classes 
comprise stars, galaxies, noise and merged objects. The catalog includes objects down to E ~ 20. 

In this study, we rejected 10 fields which had abnormal patterns of galaxies due mostly to 
spurious object detections around bright stars. We only considered the E-plate objects with 
13.5 < E < 19. This choice maximized the performance of the automatic image classifier, and 
resulted in about 530 galaxies per degree square. In this magnitude range and in our latitude 
range, the number of stellar objects outnumber galaxies by a factor of about 5. In figure |], a 
typical distribution of galaxies in this magnitude range is displayed as filled circles. The positional 
accuracy of the APM objects is better than one arcsecond. Being much smaller than our pixel 
size (64"), this uncertainty is negligible in our analysis. In order to avoid potential complications 
at POSS plate boundaries, we chose the field centers to lie further than 45' from the plate edges. 
Catalogued stellar objects found in the same fields were used as a control data set. 

It is useful to convert E-magnitudes into other magnitude systems. As was summarized by 
Koo & Kron (1992), the mean colors of galaxies is independent of magnitude for bj <; 22 and is 
close to < bj — rj > 9 ai= 1-1 i -1. Using figures 1 and 4 in Butcher & Oemler (1985), we infer that, 
for galaxies, < B-R > ga i= 1.9 ± .2, < B - V > ga i= 1-1 ± -1, and thus that < V - R > ga i= -8± .2. 
By measuring the magnitudes of stellar standards with the Automated Plate Scanner, Humphreys 
et al. (1991) have shown that, for the POSS plates, E-R~ -.011 + .148(V - R) + .058(U - i?) 4 , 
for —.2 < V — R < 1.7. Assuming that this relation also holds for the APM measurements, the 
above colors translate to < E — R > ga i= -3 ±.2, < V — E > ga i= -5 ±.3, and < B — E > ga i= 1.6 ±.3. 
Finally, figure 9 in Butcher & Oemler (1985), yields < bj — E > ga i= 1-1 ± -3. 

The mean magnitude of our galaxy sample is given by (E) = J^ max dEE^/ J^ max dE^, 
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where £ m j„ and E max are the magnitude limits of the sample, and 4§ is the mean number 
of galaxies per unit solid angle and per unit magnitude. In our sample, E m i n = 13.5 and 
Emax = 19, as stated above. A fit to the number-magnitude relation in our sample yields 
\og{dn/dE) ~ —3.44 + .30-E galaxies deg -2 (.5 mag) -1 , so that (E) ~ 17.7. However, our directly 
measured number-magnitude relation is flatter than the one measured for the southern APM 
catalog (Maddox et al. 1990), which after using the above bj — E conversion is well fitted by 
log(dn/dE) ~ —8.8 + .QE galaxies deg -2 (.5 mag) -1 , for E <^ 19. This discrepancy is probably 
due to a systematic overcorrection for saturation for galaxy magnitudes. (The slope of the stellar 
APM counts is consistent with the stellar counts of Metcalf et al. 1991). The northern and 
southern APM counts agree if we use corrected APM magnitudes defined as E' = 8.52 + .50E. In 
this case, our sample limits are E' min = 15.3 and E' max = 18.1 and the mean magnitude becomes 
(E') ~ 17.4. Given the uncertainty involved in this correction, we take our mean magnitude to be 
(E) ~ 17.5 ±.3. 

The contamination by misclassified objects in the APM galaxy catalog was estimated at 
bright magnitudes by direct visualization of 841 objects on five different POSS plates. The stellar 
contamination was observed to vary from plate to plate and as a function of magnitude, and is 
dominated by misclassified stars. The mean ratio of the number of true galaxies to the total 
number of objects classified as galaxies in the APM is {N 9 )/{N) ~ .7 ± .2 for 15 <^ E <J 16.5 
and tends to decrease with magnitude, at least in this magnitude range. Direct visualization 
for E ^ 16.5 is unreliable. Using deeper AAT plates, Maddox et al. (1990b) estimate that, 
for the southern APM galaxy catalog, (N 9 )/(N) ~ .95, for 17 < bj < 20. We thus finally set 
(N 9 )/(N) ~ .8 ± .2 as our estimate for magnitudes close to our mean (E). 



2.3. Final Sample 

From the resulting set of 67 fields, we removed overlapping fields giving preference to those 
with higher exposure times. This was done to avoid correlating the same region of the sky more 
than once. (Two fields overlapping by less than 15' were left in the sample). The final field sample 
contained Nf = 62 fields which were labeled by their Einstein sequence numbers. Their center 
positions along with the associated POSS plate numbers, exposure times, Einstein target types, 
mean X-ray intensities, and mean galaxy surface densities are listed in table |l[ Figure |l| shows one 
of the X-ray fields superimposed on the associated distribution of galaxies. 

An additional control data set was generated by scrambling the X-ray/galaxy field pairs. To 
summarize, we considered the following data sets: 

1. Hard-Galaxies: Hard-band (.81-3.5 keV) XRB (low solar contamination) crossed with 
galaxies: the real data set 

2. Hard(Solar)-Galaxies: Hard-band XRB with high solar contamination crossed with galaxies 
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3. Soft-Galaxies: Soft-band (.16-. 81 keV) diffuse emission (with a large Galactic component) 
crossed with galaxies 

4. Hard-Stars: Hard-band XRB crossed with stars 

5. Hard-Galaxies (Scrambled): As in set 1, but with scrambled field pairs. 

The first data set is the principal ("real") data set. We should not detect any correlation for the 
scrambled set (set 5), since it involves pairs of unrelated regions of the sky. As mentioned above, 
the soft band XRB is dominated by X-rays from the Milky Way. Thus, for data sets 2 and 3, 
the fraction of the X-ray counts potentially correlated with galaxies is smaller than that for the 
real set; we therefore expect the correlation to be smaller for these two sets. Stars contribute a 
non- negligible fraction of the XRB: the fraction of stars optically identified in the Einstein Medium 
Sensitivity Survey (Stocke et al. 1991) and deep ROSAT XRB surveys (see Hasinger 1996, for a 
summary) are 26% and 15%, respectively. Given the significant (mostly stellar) contamination of 
the galaxy catalog, it is important to test whether stars contribute to the CCF. This is achieved 
by the use of data set 4. After binning galaxies (and stars) into 64" pixels, the correlation analysis 
was performed in the same way for the five data sets. This provides a direct and powerful test of 
the physical and statistical significance of any resulting correlation signal. 



3. Procedure 

3.1. Correlation Function and Estimator 

Our main goal was to estimate the two-point angular cross-correlation function w xg for 
the two data sets described above. This function is a sensitive measure of correlations between 
two-dimensional data sets (see eg., Peebles 1980; Jahoda et al. 1991). To define w xg , let us 
consider two infinitesimal cells C\ and C2 with solid angles 5Q± and 5^2 with angular separation 
#12. Let 5Ni be the number of of galaxies in C\, and 5I2 be the XRB flux (in ergs cm -2 s _1 ) in 
C2. Then, w xg is defined by 

< 8N x Sh >= (n){i) [1 + w xg (9u)] 5Qi<5Q 2 , (1) 

where (n) is the mean number of galaxies per sr and (i) is the mean XRB intensity (in ergs cm -2 
s _1 sr -1 ), and the brackets denote a sky average. By taking SQi to be sufficiently small, one can 
show that equation (Q) implies that w xg (6) is the mean fractional excess in the XRB intensity at 
an angle 9 from a given galaxy. 

In practice, measurement cells have a finite size and one can not measure w xg directly. 
Instead, one computes the finite-cell estimator W xg defined as 
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where (N) and (/) are the number of galaxies and the XRB flux averaged over all available cells. 
The quantity r](9) in the numerator is defined by 

r}{9 l2 ) = {N x h)e 12 , (3) 

where the bracket denotes the product of N and I averaged over cell C\ and C 2 with a angular 
separation 0\ 2 . (Our data set, which consists of many disjoint fields, was analyzed with the 
specific averaging scheme described in §3.2.) In general, W xg depends on the cell size a and can 
thus be thought of as being a function of two variables: W xg = W xg (a, 6). If the XRB and the 
galaxy counts are uncorrelated, rj(6) = (N)(I) and thus W xg = 0. On the other hand if they are 
correlated (anticorrelated) W xg is positive (negative). 

Formally, the correlation function w xg is equal to the estimator W xg for infinitely small cell 
size, i.e. 

w xg (0) = limW xg (a,d). (4) 
Conversely, W xg is obtained from w xg as 

w xg ( a ,6 12 ) = — !L f <% 2 £ £ cm' l( m 2 w xg (e' 12 ), (5) 



Ci JC2 



where the dQ^ integrals are over two cells C\ and C 2 separated by an angle Q\ 2 with polar 
coordinates (#12, ^12) m a coordinate system whose origin is the center of C\. The d(pi 2 integral 
operates as an average over all possible relative positions of the two cells while keeping their 
separation #12 constant. We have taken the cells to be be squares with sides a. Equation (||) can 
be derived by subdividing the two cells into infinitesimal subcells and by using equation ([!]) in 
equation @. 



3.2. Implementation 

In our experiment, our minimal cells are the X-ray pixels: squares with a side length of 
a v = 64". The galaxies (and the stars) were also binned into such pixels. All invalid pixels, i.e. 
pixels under the support ribs or outside the Einstein field were ignored. We varied the cell size a 
by binning adjacent pixels together. Boundary cells, i.e. cells lying across the field edge or across 
the ribs, were rejected if less than 33% of their area were valid. The choice for this threshold was 
rather arbitrary and reflected a tradeoff between wasting pixels and having too many cells with 
large uncertainties. A moderate variation of this threshold (between 20% and 60%) did not affect 
our results significantly. We thus obtained values of W xg at zerodag for cell sizes a = na p , where 
n is an integer. We only considered one-pixel cells (a = a p ) for our non-zero lag measurements. 

For each field, we computed rj(8) (Eq. ||) as follows 

1 N c 1 N c (i,8) 
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where N c is the total number of cells in the field, Ii is the X-ray flux in d, and Nj is the number 
of galaxies in Cj. The first sum runs over all cells in the field. The second sum runs over all 
N c (i, 9) cells Cj a distance 9 away from Cj. 

This estimate is slightly different from the more standard one (see e.g. Carrera et al. 1995a, b) 
which consists in directly averaging IiNj over all pairs of cells (Ci,Cj) separated by an angle 
9. These two definitions agree in the limit of an infinite field size. For a field of finite size, our 
definition provides a consistent treatment of boundary cells. For the scrambled data set, the 
"standard" definition yielded small but significant correlation at large angle, where the finite field 
effect is most significant. These spurious correlations were absent when we used the definition of 
equation (^). We also note that, even though the roles of Ii and Nj are not formally symmetric in 
equation (^), an inversion of N and I did not noticeably affect our measurements. 

Both the X-ray and the galaxy data were subject to field-to-field variations due to differences 
in particle background contamination, solar X-ray contamination, plate sensitivity, seeing, etc. As 
a result, the fields could not be treated as a single data set. Since, in fact, we intend to focus on 
angular separations much smaller than the field size (9 <C 1°), we have instead computed W xg for 



each field separately. As we will see in §4.1, the distribution of W xg is not gaussian. However, by 
the central limit theorem, the distribution of the average of W xg over all fields is close to gaussian. 
A detailed description of the statistics of W xg is given in appendix [C]. Our final values and la 
error estimates were thus be taken be 

1 Nf 1 
W xg ± 5W xg = — J2 W xgJ ± ^^a rms (W xg ), (7) 

where W xg ,f is the value of W xg for the f'th field, and a rms is the r.m.s. standard deviation 
computed over all fields. For the scrambled data set, we constructed N scram bi e = 400 and 200 field 
pairs for the zero and non-zero lag measurements, respectively. This large number reduces the 
uncertainty in the mean and in the r.m.s. standard deviation. The final error estimate for the 
scrambled data set was normalized to Nf. 



4. Results 

4.1. Zero-Lag Results 

As we will see below, most of the correlation signal lies at 9 = 0. We therefore first consider 
the zero-lag results in great detail in order to assess the reality and uncertainty of the signal. 

Table [l| gives our measured values of W xg (a p ,0) in each field for the real data set. Figure || 
shows the distribution of W xg (a p ,0) in histogram form for the real data set. For comparison, a 
subset of the scrambled data set containing Nf fields is also shown. There is a small but significant 
positive offset in the real values which yield a correlation of W xg (a p ,0) = .045 ± .013, with a 
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significance of 3.5a. As noted above, the distribution of W xg is clearly non-gaussian for the real 
data set. This is not surprising since 77 (Eq. ||) is the product of two random variables and is 
thus expected to follow a skewed distribution. The origin of unusually high correlations for several 



outlyer fields will be discussed in detail in §4.4. Here, we simply note that, after the removal of 
the six fields with W xg > .17, the mean correlation is W xg (a p ,0) = -020 ± .008, which is now 
significant at the 2. 4a level. The scrambled set yields a correlation of W xg = —.001 ± .008 which, 
as expected, is consistent with zero. 

Figure || shows W xg (a ce u , 0) with cell sizes ranging from a p to 8a p for each of the five data 
sets. The mean values and their associated errors are listed in table Although the measurements 
at different cell sizes are not independent, their comparison gives a measure of the robustness 
of the signal. The real data set yields a positive correlation for all values of a shown. The 
solar-contaminated and Soft-Galaxies data sets yield lower correlations. As we noted above, this 
is expected since the large contribution of solar and Galactic X-rays, in each case respectively, 
results in a smaller relative contribution of the XRB to the observed intensity. 

We do not detect any significant correlation between stars and the hard-band XRB. This is 
somewhat surprising since the contribution of stars to the XRB in ROSAT deep surveys is about 
15% (Hasinger 1996), and since stars outnumber galaxies by a factor of about five in our magnitude 
range. The absence of correlation could be due to the soft nature of characteristic stellar spectra 
which mitigate their contribution in the hard band. We leave a quantitative discussion of the 
limits our analysis can place on mean stellar X-ray fluxes and on their contribution to the XRB to 
future work. The importance of this result in the context of our study of galaxies and the XRB is 
that it simplifies the interpretation of W xg . We can indeed safely assume that misclassified stars 
contaminating our galaxy sample do not contribute to the correlation signal. 

As expected, the scrambled data set yields correlation values consistent with zero for all a's. 
It is worth noting that the uncertainty for the scrambled set (5W xg (a p ,0) ~ .008) represents about 
38% of the variance (i.e. the square of the standard deviation) for the real set (5W xg (a p , 0) ~ .013). 
In general, two effects contribute to the variance of W xg : the intrinsic variation in the correlation 
signal from one region of the sky to another, and the statistical uncertainty involved in measuring 
W xg in finite fields. If we take the variance for the scrambled set to be representative of the latter 
effect, we conclude that the intrinsic field-to-field variation represents about 62% of the variance 
of W xg . One must therefore be cautious when using bootstrap error estimates since these do not 
account for the intrinsic variance. 

As a further test of the reality of the signal, we plotted W xg versus various related parameters 
in figure |4j. No obvious correlations exists between W xg , the mean number of galaxies, the mean 
XRB flux, Galactic latitude, or the X-ray exposure time. This indicates that our signal does not 
suffer from major systematic biases. 
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4.2. Non-Zero-Lag Results 

Now that we have established the existence of a significant signal at zero-lag, we can obtain 
the full correlation information by studying W xg at nonzero-lag. Figure |5] shows our measurement 
of W xg (a p , 9) as a function of 9. The results are shown both for the real and the scrambled data 
sets. The mean values of W xg and their error bars are listed in table || for each value of 9. 

As noted above, most of the signal in the real data set lies at small angles (9 < 3'). The 
scrambled set yields correlations consistent with zero, as expected. The striking feature in figure [| 
is the presence of a significant plateau for 9 ^> 4' with an amplitude of W xg (a p , 9 ^ 4') ~ .015. We 



will discuss its nature in §4.4. 



In figure || we plotted the values of W xg at 9 = 10a p against that at lla p for each field. A 
correlation (of the correlation function!) is clearly present. This shows that our measurements of 
W xg at different values of 9 are not statistically independent. This is due to the fact that the PSF 
(and perhaps the clustering scale) extends beyond one pixel. This correlation must be taken into 
account when estimating errors or when fitting to the angular dependence of W xg . In appendix [C] 
we present a treatment of this correlation using multi-variate gaussian statistics. In particular, we 
show how one can estimate the covariance matrix Vij of W xg (a p ,9i) from our data. The resulting 



values for the elements of Vij, shown in figure [13J, demonstrate that the correlation effect is most 
severe for 9 ^ 3a p and for 9 ^ 8a p . 



4.3. Comparison of Zero-Lag and Non-Zero Lag Results 

We now consider a comparison of the zero and non-zero lag signals. This is not only a check 
of the consistency of our results but also allows us to verify the reality of the correlation at large 
angle which should contribute to the zero-lag correlation for large cell sizes. 

For this comparison, we choose a phenomenological form for the underlying correlation 
function w xg which reproduces the nonzero-lag results and check that it is consistent with the 
zero-lag results. We take w xg {9) to be the sum of a gaussian and a constant, i.e. 

w xg (9) = a + ai e-^ d / e ^, (8) 

where oo, a\ and 9q are constants. The estimator W xg can then be computed by introducing 
equation (||) in equation (||). We obtain 

W xg (a,9) = a + a 1 A(a,9 ,9) (9) 



where A is the integral defined in equation (Bl). A semi-analytical expression for A is given in 
appendix BA and its angular dependence is plotted in figure 12 for a = a p and 9q = 0'.5. 

Our best fit parameters for W xg (a p ,9) at nonzero-lag are ao = .013 and a\ = .087. For 
convenience, we fixed 9$ = a ps f = 0'.5. The result of the fit is shown in figure || as the dot-dashed 
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curve. The resulting cell size dependence of W xg at zero-lag is shown in figure ||, as a similar curve. 
The agreement is acceptable even though the measured values of W xg are slightly higher than the 
one inferred from the non-zero lag results. This is probably due to the ambiguity in the treatment 
of boundary cells (see discussion in § |3.2[ ). (In fact, a variation of the pixel inclusion threshold can 
slightly improve the agreement). Note also that the zero-lag measurements confirms the presence 
of a significant correlation all the way up to a ~ 8', which corresponds to 9 ~ 4'. 



4.4. Correlation at Large Angles 

The plateau in W xg at large angle is rather surprising. Its reality is confirmed both by the 
scrambled results and by the fact that it appears consistently both in the zero and non-zero lag 
measurements (at least for 9 ~ 4'). 

In previous studies (Lahav et al. 1993; see also the other cross-correlation studies discussed 
in §|l]), the correlation signal was taken to arise from two terms in the correlation function: the 
"Poisson" term and the "clustering" term. The Poisson term corresponds to the case of the 
galaxies themselves emitting X-rays, and the clustering term corresponds to the spatial clustering 
of galaxies with discrete X-ray sources. As we will see in detail in §[|, these two terms produce, 
in our context, correlation signals with specific angular dependence given by the two functions A 
and 3 displayed in figure [l^. For both terms, W xg falls off for 9 ^ 4' and thus can not produce 
the plateau observed in our measurement. This can be seen more quantitatively in figure [| where 
the dotted curve correspond to the best fit to the Poisson + clustering terms (see §^3). The fit is 
poor for 9^4'. 

One possibility is that the plateau is due to the hot X-ray gas in clusters of galaxies. Diffuse 
emission from the clusters in our fields is mostly, but not completely, removed by our point 
source excision algorithm and can thus contribute signal to the very sensitive CCF. An example 
is provided by the field displayed in figure [7]. This field has the highest zero-lag value of W xg in 
our data set (i.e. it is the right-most field in figure §). A cluster of galaxies appears as an excess 
in residual X-rays and in the surface density of galaxies in the upper left-hand corner of the field. 
Bright nearby clusters have typical X-ray sizes of about 10', a scale which is consistent with the 
scale of the plateau. Note that the field displayed in figure ^ is extreme; most fields in our data 
set, including those with cluster targets, such as the one in figure ffl, do not show any obvious 
cluster emission. A search of the NED database (Helou et al. 1991), indeed reveals that only three 
of the cluster-target fields have an X-ray detected cluster within a radius of 20' of the field center. 

To test this hypothesis we use the fact that 29 out of our 62 fields had clusters of galaxies 
as their Einstein targets. These are flagged with a target type of 3 in table |l[ (Fields with 
superclusters as targets often contain clusters of galaxies and were thus also included). In figure || 
we show the value of W xg at 9 = vs. that at 9 = 7a p with different symbols for fields with 
cluster and with non-cluster targets. The value of W xg (a p , 7a p ) is taken to be representative of 
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the plateau for each field and is listed in table ||. As can be seen in this figure, the zero-lag value 
has little relation to the presence of a cluster, except for a few outlyers. This is in contrast with 
the non-zero lag values which are significantly larger for cluster target fields. The values of W xg as 
a function of 9 for each of the cluster and non-cluster data sets are listed in table ||. At zero-lag, 
the non-cluster fields yield a mean correlation of W xg (a p ,0) = .029 ± .013, a significance of 2.2a. 
In figure g the mean W xg is shown separately for cluster and non-cluster fields. The value for all 
fields combined is redisplayed for comparison. In spite of the somewhat degraded statistics, it is 
clear that the plateau is virtually absent for non-cluster target fields, whereas it is enhanced for 
cluster target fields. This suggests that the plateau is indeed caused by the residual diffuse X-ray 
emission in clusters of galaxies. 



5. Interpretation of the Correlation Function 

Now that we have established the existence of a positive correlation signal for our data 
set, we need to interpret its angular dependence and thus infer information about galaxies and 
the XRB. For this purpose, Treyer & Lahav (1995) have developed a detailed formalism which 
involves a projection of the 3-dimensional CCF into 2-dimensions. Here, we consider a similar but 
fully 2-dimensional approach. This has the advantage of being simple and of involving few model 
parameters. In addition, we include the effect of the PSF of the X-ray imaging detector. The 



formalism presented in §5.1 is general and can be adapted to any measurement of the CCF between 
an imaged diffuse background and a population of discrete sources. After applying this formalism 
to our specific case in jj|5.2|, we compare its prediction with our non-zero lag measurements in §|5.3| 



5.1. Formalism 

In general, the total X-ray flux Ij = Ij measured in a cell Cj can be decomposed into 

if = ij + q + if, (io) 

where the terms, correspond to the instrumental background (i) (due to cosmic ray interactions, 
calibration source leakage, and solar-scattered X-rays), the emission from point X-ray sources (p), 
and the diffuse component of the XRB (d), from left to right respectively. A source is considered 
discrete if its solid angle is negligible compared to the PSF and to the cell size. Thus, in our 
situation, the diffuse component of the XRB includes the emission from the intracluster gas in 
nearby clusters of galaxies. We can also decompose the galaxy catalog object counts Nj = N° in 
Cj as 

Nj = N 9 j + N s j (11) 

where Nj is the number of galaxies (g) in Cj and Nj is the number of spurious objects (s) 
misclassified as galaxies. 
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introducing equations (10) and (11) in equation (|3|), we obtain 

r] = rjis+ ilps + Vds + Vig + Vpg + Vdg ( 12 ) 

where ij a },(6) = (IfNng i:j= g. As before, the brackets refer to an average over cells separated by an 
angle 8. 

Since the galaxy catalog and the X-ray maps were constructed independently, P and N° are 
not correlated. Thus, r] is + rj ig = (P)(N 8 ) + (P)(N 9 ). 



As we noted in §4.1, stars are not significantly correlated with the XRB intensity, at least 
in our experimental conditions. Other misclassified objects such as plate flaws are obviously not 
correlated with the XRB. We thus obtain r] ps + r] ds = (I P )(N S ) + (I d )(N s ). 



We have seen in § [4.4| that the diffuse component of the XRB is likely to affect the CCF. An 
analysis of the diffuse term would involve a modelling of clusters of galaxies and of their redshift 
distribution. This is beyond the scope of the present study. We therefore leave rjd g unexpanded. 
In §5.3, we will show how we can nevertheless proceed using the subset our measurements for 
non-cluster target fields. 

We are left with r] pg {8) to evaluate. For this purpose, we define the PSF, ijj, so that a point 
source with flux S p , placed at P , produces an observed X-ray intensity at 9 given by 

i{0 o ) = S P V(M- ( 13 ) 

We assume that ip{6) is azimuthally symmetric so that it only depends on 8 = \6\. This definition 
implies that ip is normalized so that 



dn a ii>(9 a ) = i, (14) 



where the integral is over the whole sky. 



The derivation of r] pg using this definition of ifi is given in detail in appendix The result 
consists of the two standard terms (see eg. Treyer & Lahav, 1995) modified to include the effect of 
the PSF: the "Poisson" term (Eq. |A5| ) and the "clustering" term (eq. [|A10[| ). Introducing these 
results into equations (|i~2|) and @, we finally obtain 

W xg = W"P 9 oisson + W^ ustering + W*f use , (15) 

where 



w ™ ( 12 ) = A_L __ / #12 / dni / dn' 2 W 12 ) (16) 



{19) 1 
(J) (AT) 27ra 2 Jo ^ Jc^'" 1 Jc 2 
and 

^clustenng^^ = J_ Fjfa f ^ f ^ f dO^^^O . (17) 



(I)(N) W7o ^ jcr 'Jci 2 Js s ~ p9K 

In the two equations above, the dVt'i and d£l' 2 integrals are over two cells C\ and C2 separated by 
012- The dtt' s integral is over the whole sky (S). The brackets denote an average over all available 
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cells. The function w pg (6) is the effective CCF between the positions of point X-ray sources and 
galaxies. It is defined in equation ( |A9[ ) as the X-ray flux-weighted average of the flux-dependent 
CCF w P g. That such an averaging is necessary is not surprising since, after all, the X-ray point 
sources with fluxes just below the excision threshold dominate the anisotropy of the observed 
XRB and are thus the one most susceptible to producing a correlation with the galaxies. As noted 
above, the term W* g Susc = (r] dg - {I d ) (N 9 ))/((/) (JV» is left unexpanded. 

As can be inferred from inspecting their functional form, the three terms above correspond to 
three distinct effects responsible for a correlation between the XRB and galaxies: 

1. Poisson: the galaxies emit X-rays. 

2. Clustering: the galaxies are spatially clustered with discrete X-ray sources. 

3. Diffuse: the galaxies are spatially clustered with diffuse X-ray emission. 



5.2. Application 

We can now apply the above formalism to our specific case. For the Einstein-IPC, the PSF is 
gaussian so that, after normalization (Eq. flli}), 



= ^^-e-^Aw) 2 (18) 



To proceed, we need to know the angular dependence of W pg which describes the clustering of 
galaxies with point X-ray sources. Information about the clustering of sources can be deduced from 
the measurement of their Auto-Correlation Function (ACF). The angular ACF w gg of galaxies 
in the southern APM catalog is well described by a power law of the form w gg (6) = (6 l®gg) ^ 99 
in the range .01° < 9 < 1° (Maddox et al. 1990a). The slope was measured to be (3 gg = .668 
and is magnitude- independent. On the other hand, the correlation amplitude 8 gg does depend on 
magnitude and ranges from about 10" to 151" for Bj ~ 20 and 17.5, respectively. Such a scaling 
of the ACF with magnitude is expected from the projection the spatial ACF into 2-dimensions 
(Groth & Peebles, 1997; Phillips et al. 1978). Vikhlinin & Forman (1995) studied the clustering 
of ROSAT X-ray sources with a .5-2 keV flux lower limit of about 10~ ergs cm~ 2 s^ 1 , close to 
our flux excision limit S exc . They found that the ACF of ROSAT point X-ray sources, w pp , also 
follows a a power law with a slope of (3 PP = .7 ± .3 and a normalization of 6 PP = 4" for 9 between 
30" and 1000". 

The two slopes (5 gg and f3 pp are both consistent with the value ~ .7. It is thus reasonable 
to assume that the magnitude-dependent CCF, w pg , defined in equation (|A7[ ) is also described 
by a power law of the form: w pg (S p , S 9 ;6) = (8/9 pg (S p , S 9 ))~@, where S p and S 9 are the X-ray 
and optical fluxes of the point X-ray sources and galaxies, respectively, and (3 ~ .7. Since the flux 



-16- 



dependence of w pg is limited to its normalization 9 pg (S p , S 9 ), we can factor the 9 dependence of 
w pg out of the integrals in equation (A9) and obtain 

Wpg{6) = ' (19) 

where 9 pg is left as a free parameter to be determined. 



With the specific functional forms of equations (|lq) and (19), equation (|iq ) becomes 

W%°™™(a,e) = ^ 2 Ma,a psf ,9), (20) 

while equation ( |T7| ) becomes 

W^ ste ^( a ,9) = V'W 9 } Z(a,a psf ,f3,9), (21) 



'P9j 



where A and H are the integrals defined in equations (|B1| ) and (|B8|), respectively. In Appendix [b] 
we describe our semi-analytic evaluation of A and our numerical evaluation of E. The functions 
are both displayed in figure [l2|. Two useful analytical fits to A and S for the parameters relevant 
to the present study are given in equations (B7) and (B9), and are also displayed on this figure. 
(The analytical fit for A is indistinguishable from A itself on the figure.) 



5.3. Fit Results 

Our remaining task is to compare our theoretical expectations with our measurement of the 
CCF. The three terms contributing to W xg have different angular dependences and can thus, at 
least in principle, be isolated by a fitting to measurements of W x „ at different values of 9. 



As we noted in £4.4, the Poisson and clustering terms can not produce the plateau observed 
for the combined cluster plus non-cluster data set. The best fit parameters to a combination of 
the poisson term and the clustering term (see Eq. p3] below) are a po isson = and a c i ustcr ing — -05. 
(The two parameters are constrained to be positive). The best fit is shown in figure [5] as the 
dotted curve. As we already noted, the fit is poor for 0^4'. 

Given the disappearance of the plateau when only non-cluster target fields are considered (see 



§4.4), it is likely that the plateau is due to residual diffuse X-ray emission from clusters of galaxies 
which contribute to W xg through the diffuse term. Although we do not pursue this here, it would 
be worthwhile to model such diffuse emission and to derive the expected angular dependence of 

^/diffuse 
xg 

Let us now focus on the non-cluster target fields for which the plateau is virtually absent. 
If we assume that the diffuse term is negligible for the non-cluster data set, we are left to fit 
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W X g(ct,6) to the following function 

W xg (a, 9) = a poisson A(a, cr ps f,0) + M, (22) 

where a po isson) and a c i uster i ng are free parameters. So as to reduce the number of free parameters, 



we fix (3 = .7 in accordance with the discussion in §5.2. 



We have noticed in §£hj, that the measurements of W xg at different values of 9 are not 



statistically independent. In appendix^, we describe our procedure for the fit (Eq. |C5|) assuming 



that the set of values of W xg follow a multi-variate gaussian distribution (Eq. []C1| ) . The covariance 



matrix (Eq. [p2fl ) is estimated from the data (Eq. \C4\) and is displayed in figure |T 



The result of the fit to our measurement of W xg (a p ,9) for < 9 < 8a p are displayed in 
figure 10. This figure shows the constant probability contours of x 2 (Eq. |C7f| ) with 9 — 2 = 7 



degrees of freedom. We have used the analytical fits for A and S given in equations ( [B7D and (B9) 



respectively. The best fit (Eq. |C8| ) occurs for a po isson = -04 ± .02 and a c i us t e ring = -005 ± .004. 
The errors correspond to the la estimates for each parameter taken separately (Eq. [ |C10| ) and 



are also displayed on figure 10. The resulting best fit is displayed as the solid line in figure 



The value of x 2 f° r the best-fit parameters is 5.1. The fit is thus acceptable. It worth noting 
that a fit which does not takes the covariance of W xg into account (i.e. one for which the covariance 



matrix V, in Eq. [02], is diagonal) yields x — 2-1 which correspond to a fit which is "too good". 



The covariance of W xg must therefore be taken into account in order not to underestimate the 
measurement errors. 



By comparing equations (20) and (^) with equation (^), we can deduce the following 
physical quantities 



(19) = ^|£/(iV)(/)a poisson (23) 
or 



and 



(I)(N) 



i 



U P9 ~ a ^ /Jp\/JV9) ) "clustering" ( 24 ) 

Our measurements for the non-cluster data set yield (N) = .149 db .002 objects per cell 
and (J) = (3.11 ± .01) x 10~ 4 cts s _1 per cell. From our estimates in section 0, we have 
(N 9 ) /(N) = .8 ± .2 and (I XRB )/(I) = .6 ± .1. If we assume that the diffuse component of the XRB 
is negligible, we can set (I p ) ~ (I XRB ). (Alternatively, any difference between (I p ) and (I XRB ) 
can be absorbed in the definition of 9 pg ). With these numerical values and setting a = a p , we 
obtain (i 9 ) = (2.2 ± 1.1) x 10~ 6 cts s -1 arcmin -2 and 9 pg = .10 ± .16 arcsec. The errors correspond 
to la. The implication of these final values are discussed in the next section. 



6. Galaxies and the X-Ray Background 



-18- 



6.1. Clustering 

Let us first consider our measurement for the galaxy /X-ray-sources clustering amplitude 9 pg . 
Our value for 9 pg is consistent with at the la level. Note that negative values of 9 pg are of course 
unphysical and are thus excluded. Therefore, we do not detect any significant clustering between 
discrete X-ray sources and galaxies. Our value for 9 pg can thus be interpreted as an upper limit. 



As we noted in section 5.2 , the amplitudes for the auto-correlation of galaxies (Maddox 
et al. 1990a) and of ROSAT X-ray sources (Vikhlinin & Forman 1995) are 10" < gg < 151" 
(for 20.5 > bj > 17.5), and 9 PP ~ 4", respectively. Our value 9 pg = .10 ± .16 arcsec is thus 
considerably smaller than these two amplitudes. In our analysis, we are in effect cross-correlating 
two population of sources with different mean X-ray flux: the X-ray sources dominated by the 
sources just below the excision flux limit S exc (i.e. mostly AGN at z > .2), and the galaxies 
(at z .2), which, as will will see below, have a much smaller mean X-ray flux. Thus, it is 
not surprising that our cross-correlation clustering amplitude be smaller than the two associated 
auto-correlation clustering amplitudes. 



6.2. X-ray Emission from Galaxies 

The statistical significance of our detection for the Poisson term is also low: our measurement 
of (i 9 ) is only 2a away from 0. This can easily be improved by performing the analysis with a larger 
number of fields. We can nevertheless draw several useful conclusions from our measurement. 

Using our estimation of (I XRBs j (see section ||), we deduce that the fraction of the XRB 
contributed by APM galaxies with 13.5 < E < 19 is (I 9 )/{I XRB ) = .015 ± .008 in the hard (.81-3.5 
keV) band. This fraction should not be mistaken for the total contribution of APM-like galaxies 
distributed over all space. 

The mean X-ray flux of the catalogued galaxies is easily computed to be 
(S 9 ) = (i 9 )/(n 9 ) = (2.1 ± 1.1) x 10 -5 cts s _1 in the hard band. The X-ray spectra of galaxies can 
described by a thermal bremsstrahlung model with kT ~ 1 — 2 keV and kT ^ 5 keV, for ellipticals 
and spirals, respectively (Kim et al. 1992). If we adopt such a model with kT = 3 db 2 keV and 
N H ~ 3 x 10 20 cm" 2 , the above flux corresponds to (S 9 (.8l - 3.5keV)) = (8.1 ± 4.7) x 10~ 16 
ergs cm -2 s . In the more conventional Einstein and ROSAT bands, this corresponds to 
(S 9 (.2 - 4.0keV)) = (9.1 ± 5.0) x 10~ 16 and (S®(.5 - 2.0keV)) = (5.4 ± 3.0) x 10~ 16 ergs cm" 2 
s , respectively. We note that this flux is three orders of magnitude below the flux limit of the 
Einstein Medium Sensitivity survey (~ 10~ 13 ergs s _1 cm -2 in the .2-4 keV band; see Gioia et al. 
1990) and more than one order of magnitude below those for deep ROSAT surveys with mostly 
complete optical identifications (~ 10~ 14 ergs s -1 cm~ 2 in the .5-2 keV band; see Hasinger 1996). 
This is of course a consequence of the statistical nature of our measurement. 



It is noteworthy that (S 9 ) is about two orders of magnitude below the discrete source 
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detection threshold in our fields ((S exc (.81 — 3.5keV)) ~ 4.8 x 10 -14 ergs cm -2 s _1 ). Even though 
this seems surprising at first, we will see below that (S 9 ) is consistent with the mean X-ray flux 
expected for galaxies in our sample, i.e. for galaxies with (E) ~ 17.5. On the other hand, galaxies 
with X-ray fluxes close to S exc are expected to have a mean magnitude of (E) ~ 13.5. They are 
thus, on average, much brighter in the optical than the galaxies we considered. Another way to 
understand the low value of (S 9 ), is to consider the surface density of X-ray sources. According to 
the X-ray logN-logS relation measured with ROSAT (Hasinger et al. 1993), the number density 
of X-ray sources with fluxes above (S exc ) is n(> (S exc )) ~ 7 deg -2 , much lower than our galaxy 
surface density ((n 9 ) ~ 530 deg -2 ). On the other hand, an extrapolation of this logN-logS relation 
yields n(> (S 9 )) ~ 10 3 deg~ 2 , close to (n 9 ). 

The X-ray emission of "normal" galaxies has been extensively studied in the past (see 
Fabbiano 1989 for a review). For spiral (S) and irregular (Irr) galaxies, the X-ray emission is 
dominated by binary X-ray sources and SNR. For elliptical (E) and SO galaxies, it is likely to be 
composed of a baseline component of X-ray binaries and a hot gaseous component. Fabbiano et 
al. (1992) have constructed a catalog of Einstein X-ray measurements for more than 450 galaxies. 
The catalog is mostly optically selected and is thus appropriate for a comparison with our results. 
We selected the 222 and 132 galaxies with Einstein-IP C fluxes classified as "normal" S+Irr and 



E+SO galaxies. In figure 11, we show the resulting S(. 2-4.0 keV) X-ray fluxes vs. the optical 
B-magnitudes for S+Irr and for E+SO galaxies separately. Flux upper limits are indicated by 
stars. Early-type galaxies tend to have higher X-ray fluxes than late-type galaxies. This is likely 
due to the contribution of the hot gaseous component which is absent in most late-type galaxies. 

For S+Irr galaxies, a correlation between log(S) and B is significant at the 99.99% confidence 
level (ignoring upper limits). The correlation for E+SO galaxies is only significant at the 91% 
level. These confidence levels were computed using the linear correlation coefficient for detected 
galaxies. If, as is customary, we assume a constant X-ray-to-optical flux ratio we obtain the 
relation log S = —\B — c, where c is a constant. Ignoring galaxies with upper limits for S, we find 
c = —7.9 ± .6 and c = —7.3 ± .6, for the late and early type galaxies, respectively. The errors 
correspond to la. (The above late-type relation is only slightly different from the one proposed 
by Fabbiano et al. which corresponds to c = —8 ± 1). These relations are shown as the dashed 



lines in figure 11|. Note that the early-type relation must be taken with caution since no strong 
correlation is found for this sample. In general, the fraction of late-type to early-type galaxies 
must be known to predict the mean flux of our mixed APM sample. In the absence of a more 
significant correlation between X-ray and optical flux for the early-type sample, a comparison with 
each sample separately is sufficient. 



Our measurement is shown as the cross in figure [Tl|. The conversion of §2.2 was used to 



transform (E) into (B) = 19.1 ± .4. As is apparent on the figure, our results agree within la with 
the X-ray-to-optical flux relations inferred from the catalog of Fabbiano et al., simultaneously for 
both the early and late type samples. We conclude that our correlation results are consistent with 
an extrapolation of the known X-ray-to-optical flux ratio of normal galaxies. Our result therefore 
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confirm earlier estimates of the contribution of galaxies to the XRB obtained by integrating the 
optical number-magnitude relation of galaxies combined with the assumed constant X-ray to 
optical flux ratio. By applying this method to earlier but similar Einstein measurements of the 
X-ray-to-optical flux ratios (Fabbiano & Trinchieri, 1985; Trinchieri & Fabbiano, 1985), Giacconi 
& Zamorani (1987) concluded that galaxies contribute about 13% of the 2 keV XRB. Our results 
do not exclude the possibility that a change in the X-ray-to-optical flux ratio for galaxies with 
E ;> 19 results in a larger contribution to the XRB. 

7. Conclusions 

Our measurement of the cross-correlation function estimator W xg between 62 Einstein-IPC 
fields (.81-3.5 keV) and 13.5 < E < 19 Northern APM galaxies therefore leads to the following 
conclusions: 

At zero-lag, we detect a 3.5c correlation signal with an amplitude of W xg (a p , 0) = .045 ± .013. 
This signal passes our series of control tests. Intrinsic field-to-field variations contribute about 62% 
of the variance of W xg . We therefore warn against the sole use of bootstrapping error estimates 
which do not include this source of uncertainty. 

At non-zero lag, the angular dependence of W xg (a p , 9) has two main features: the main signal 
for 9 <; 4' and, for 9 ^ 4', an almost flat plateau with an amplitude of about W xg (a p , 9 ^ 5') ~ .015. 
When fields whose Einstein targets are cluster of galaxies are removed from our sample, the 
plateau virtually disappears, and the zero-lag amplitude becomes W xg (a p ,0) = .029 ± .013. The 
agreement between zero-lag and non-zero-lag results were shown to be acceptable. 

We developed, from first principles, a formalism to interpret the angular dependence of W xg 
taking into account the PSF of the X-ray imaging detector. We find that the correlation signal 
can be produced by three distinct effects: 1. the X-ray emission from the galaxies themselves, 
2. the clustering of galaxies with discrete X-ray sources, and 3. the clustering of galaxies with 
diffuse X-ray emission. These three contributions correspond to the Poisson, clustering and diffuse 
terms in the expansion of W xg . These terms have different angular dependences, and can thus, in 
principle be isolated with our formalism. 

Even though the diffuse component probably contributes little to the total intensity of the 
XRB, it can have a significant effect on W xg through the diffuse term. It is likely that the 
diffuse X-ray emission from clusters of galaxies only partially excised by our point-source finding 
algorithm produces the oberved plateau in W xg . 

After assuming that the diffuse term is negligible for our non-cluster target subsample and 
that the CCF between galaxies and X-ray point sources has the form w pg (9) = (9/9 pg )~' 7 , we could 
fit W xg (9) with the Poisson+clustering term. For the fit, we took proper care of the statistical 
correlation of measurements of W xg at different values of 9 which were shown not to be statistically 
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independent. We find that the galaxy X-ray sources correlation amplitude 9 pg is .10 ± .16 arcsec. 
As expected, this is smaller that the amplitude for the autocorrelation of galaxies and X-ray 
sources separately. Our value for 6 pg is consistent with zero at the la level and can thus be taken 
as an upper limit; we therefore do not detect a significant clustering between galaxies and point 
X-ray sources. 

From the consideration of the fields with non-cluster targets, the mean X-ray intensity 
produced by galaxies in the .81-3.5 keV band is (i 9 ) = (2.2 ± 1.1) x 1CT 6 cts arcmin" 2 s" 1 
which corresponds to 1.5 ± .8% of the XRB in this band. The mean X-ray flux of galaxies with 
(E) = 17.5 ± .3 is then (S 9 {.81 - 3.5keV)) = (8.1 ± 4.7) x 1(T 16 ergs s" 1 cm" 2 . This is well below 
the flux limits of both Einstein and ROSAT surveys with mostly complete optical identifications. 
A comparison with the X-ray catalog of galaxies of Fabbiano et al. (1992) shows that this flux 
agrees within la with that expected for both early and late-type galaxies assuming a constant 
X-ray-to-optical flux ratio. With such an assumption, Giacconi & Zamorani (1987) concluded that 
normal galaxies contribute about 13% of the 2 keV XRB. Our results, however, do not rule out 
the possibility of a larger contribution from galaxies with E ^ 19. 

From an observational standpoint, it would be worthwhile to improve the statistics by 
increasing the number of fields. It would also be useful to apply the technique with high statistics 
to deeper galaxy samples and to other X-ray bands. From a theoretical point of view, the diffuse 
term and, in particular, the effect of diffuse emission from galaxy clusters deserve more attention. 
Investigations of the non-gaussian distribution of W xg would allow a better treatment of errors 
and would perhaps allow the extraction of more information from the correlation signal. 

It is important to note that our formalism does not require any input parameters. (With 
better statistics, the exponent for w pg , which was fixed to (3 = .7 in this study, can be left as a free 
parameter). Results for (S 9 ) and 8 pg can therefore be considered as measurements rather than 
as a fit to a model. This cross-correlation technique is thus a powerful way to measure the mean 
flux and clustering amplitude of sources too faint to be resolved. In addition, if our conclusions 
regarding the origin of the plateau are confirmed, this technique can perhaps be used to probe 
faint diffuse X-ray emission which would otherwise be undetectable. 
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Caltech, under contract with the National Aeronautics and Space Administration. This work was 
supported by the National Aeronautics and Space Administration under the Long Term Space 
Astrophysics Research Program grant NAGW2507. This paper is Contribution Number 602 of the 
Columbia Astrophysics Laboratory. 
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A. Computation of r) pg 

In this section we compute r] pg , the main ingredient in the interpretation of W xg , taking into 
account the effect of the PSF of the X-ray instrument. This quantity was defined in Section [Bj] as 

n Pg (e) = (Jfiv/) ei . =e , (Ai) 

where if is the flux in Cj from point X-ray sources, N? is the number of galaxies in Cj, and the 
average is over cells with separation 6. 

Using the definition of the PSF (eq. [jOJ), we can write 

if=[ dQi [ dn s P(e a )^(0 si ), (A2) 

Jd Js 

where i p (0 s ) is the intrinsic intensity of the XRB due to point sources at 6 S (i.e. the intensity 
which would be measured with a delta- function PSF). The second integral is over the whole sky. 

We proceed in the usual way (Treyer & Lahav 1995), by partitioning the sky 5 and the cell 
Cj into infinitesimal subcells labeled by k and I, respectively. Then, r\ vg becomes 



Vpg 



[ dQ, £ 5P k m k ) ) ( E SN l ) ) (A3) 
Ci fees / Wc, / / 



where 51^. is the intrinsic X-ray flux from point sources in subcell k, anf is the number of 
galaxies in the subcell 1. 

Upon isolating the instances where the k and / subcells coincide, we can separate the above 
equation into two terms, namely 



Vpg 



' J c dn, ( £ * ijav* + E siIsnA v>(*«)\ (A4) 



k^l J 

Following common us£ig6 (s66 6g. TVcycr Sz La-liav, 1995), we denote these two terms by ?7poisson 
and Clustering, respectively. 



A.l. Poisson Term 

If the k-subcells are sufficiently small, they will contain either or 1 galaxies. As a result, 
the quantity inside the bracket in ?7poisson is equal either to or to the X-ray flux of the galaxy, 
respectively. Upon replacing the k-summation by an integral, we thus obtain 

%,oisson(%) = (I 9 )^ dfa J c dn[ J c cKfyKfly). (A5) 

where {I 9 ) is the mean X-ray flux emitted by the catalogued galaxies into a cell of side length a. 
As usual, the dfaj integral operates as an azimuthal average over the relative positions of the two 
cells. 
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A.2. Clustering Term 

The treatment of the clustering term is more complicated. This is due to the fact that this 
term involves, as its name indicates, the clustering between X-ray sources and galaxies, which, in 
general, is flux-dependent. Let us thus decompose the product SI^SNf in ?7 c i us tering (eq. |A4fl ) as, 

( A6) 

where S p is the X-ray flux of the X-ray point sources, and S 9 is the X-ray (or optical) flux of 
the galaxies. The quantities d(5N%)/dS p and d(5Nf) / dS p are, respectively, the number of X-ray 
sources and galaxies in the k and 1-subcell per unit flux. To simplify the notation, the integration 
limits relevant for both data sets, although present, have not been displayed explicitly. 

For sufficiently small subcells, the product in the bracket is equal to the joint probability 
8P^f(S p , S 9 ) of finding an X-ray source in the k-subcell with flux S p and a galaxy in the 1-subcell 
with flux S 9 . This probability can be parametrized in terms of the CCF w pg between X-ray 
sources and galaxies in the usual fashion (see eg. Eq. Q) as 

8P™(S p ,S 9 ) = [l + tito(^,5»;*u)]*n*«l, (A7) 

where 5Vt m is the solid angle of subcell m, and n refers to the appropriate number of sources per 
unit solid angle. The flux dependence of w pg is explicitly displayed. 

Equation (A6) can thus be conveniently written as 

(si^SNf) = {i?)(n 9 )[l +W pg {9 k i)]Sn k 5n l , (A8) 

where (i p ) = J dS p S p (dn p / dS p ) is the mean X-ray intensity due to point sources, and the effective 
CCF Wpg is defined by 

s m^)S dS 'l dsa Q Q^^^)- < A9 > 

With this in hand, and by turning summations into integrals in the second term of 
equation (|A4|), we finally obtain 

1 r 2jT 



%lu S tering(012) = {I P )(N 9 ) 



l+ 2^U #12 l Cl ^ ic ^ is dn 's Wp9 ^)^( 2s) 



(AW) 



B. Useful Functions 

B.l. Semi-analytic Calculation of A 

In this appendix, we evaluate the integral A(a, <r, #12). This expression is handy for our 
phenomenological fit of W xg (§ [4,3| ) and appears in the Poisson term of W xg when the effect of 
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the PSF is taken into account (§ |5.2| ). Let us consider two 2-dimensional square cells, C\ and C 2 , 
with mutually parallel sides of length a. Denote their separation by the vector #12 whose polar 
coordinates are (#12, ^12) i n a coordinate system with origin at the center of C\. Then, A is defined 
by 

1 r 2n 

A(a,a,6 12 ) = — #i 2 A(a, a, 6*12,^12) (Bl) 
2ir Jo 

where the d(f)\ 2 integral operates an azimuthal average of A itself defined by 

k(a,a,6 12 Ai2) = —J d% [ d 2 6' 2 e~^l^\ (B2) 
a JCi Jc 2 

where the d 2 9 integrals run over the two cells, and 8' 12 = — 0L\- The a -4 factor was included 
to make A dimensionless. 

In general, the three length scales involved in A (namely a, 9\ 2 and a) are comparable 
and, thus, one can not apply any simplifying approximations. The integral can nevertheless be 
performed analytically. For this purpose, let us choose the axes of the above coordinate system to 
be parallel to the cell sides. One can then conveniently separate A as follows 

A = ^jAi 2 (a, a, x 2 )A 12 (a, a, y 2 ), (B3) 
or 

where (x 2 ,y 2 ) are the coordinates of the center of C 2 and, 

pa/2 rh+a/2 . 

A(a,a,t 2 )= dtx dt 2 e-^ t2 ~ ll) /<J . (B4) 

J-a/2 Jt 2 -a/2 

Note that, in this notation, x 2 = 0i 2 cos(4>i 2 ), and y 2 = 6i 2 sm((pi 2 ). We operate a change of 
coordinate defined by: t\ = t 2 — t\, t 2 = t\ + 1 2 . The integration region in Eq. B4 turns from a 
square to a diamond which can be subdivided into four regions. The symmetry about the t 2 = i 2 
axis allows integration inside the upper two regions only. Thus, A becomes 

1 / rfa rn+a rh+a r~ri+2t2+a \ . 

- x 2 / dn dr 2 + dn dr 2 e~3^/ ff ) . (B5) 

1 \Jt 2 -a Jt 2 Jt 2 hi ) 

The \ factor comes from the Jacobian of the coordinate transformation. 



After performing the two integrals in Eq. B5 and after some algebra, we finally get 
A(a,a,t) = a 2 



e 2 



FX r ~ f t \ , . „ ft — a, 

+ W-a[-2terf — + (t - a) erf 



2 \2a J v V 2a 

+ c 
2a 



+ (t + a)evf[ t -^-)} (B6) 



where erf is the error function. 
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To evaluate A we need to average the above result over <f>i2- We do so by numerically 



averaging A over 20 values of 0i2 between and 7r/4. Figure 12 shows the resulting values for 
A(9) as a function of 6 for the parameters relevant to our study, i.e. for a = a p = 64" and 
a = (Tpgf = .5'. For comparison with other relevant functions, A(0) was normalized to 1. The 
actual normalization is A(0) ~ .551. 

It is convenient to fit the following analytical form to A(6) 

A(0) ~ Aie"^^ 2 ) 2 , (B7) 

where Ai and A2 are coefficients whose best-fit values are .552 and .677 arcmin, respectively. The 
fit results in a mean residual RMS of .0052 and is indistinguishable from the solid line in figure O. 



B.2. Numerical Evaluation of 3 

The other useful function, 3, appears in the clustering term of W xg when the effect of the 



PSF is included (§p\2|). Let us consider again the two cells, C\ and C2, described in appendix |B.l 
The desired function 5 is defined as 

E(a,a,p,e l2 ) = 2 \ 2 / 2 "# 12 / dn[ f dn' 2 f dn s e -K^/-) 2 (B 8) 

where 6' is = \0'- — 9' |, i = 1, 2. As usual, the d(j)\2 integral performs an azimuthal average. The 
dQ[ dO,' 2 integrals are over C\ and C2, respectively. The dO, s integral is over the whole sky. The 
constant in front of the integral was included to make 3 dimensionless. 



Given the cumbersome nature of equation (B8), we approximated 3 by Monte-Carlo (MC) 



integration for the parameter values relevant to our analysis, i.e. for a = a p = 64", a = a ps f = .5', 
and (3 = .7. The desired accuracy was set to .04, and the all-sky integral was approximated by a 
square cell of side 25a p centered on the midpoint between C\ and C2. 

The result of the MC integration is shown in figure [l^ for values of 9 ranging from to 8'. As 
before, 3 has been normalized to 1 at 6 = 0. Its actual normalization is 3(0) ~ 1.480. 

It is here also convenient to fit 3(0) with an analytical function to facilitate modelling of the 
observed correlation function. We choose the form 

3.0) - [^-^j CM) 
where £0 is the only free parameter. The best fit occurs for £0 — -573 ± .008 arcmin and yields a 



reduced x 2 °f -32 for 28 degrees of freedom. The result of the fit is also shown in figure 12. The 
power law (8/a)~^ (after being renormalized) is also plotted and can be seen to agree well with 5 
for 8 ^ 3/. This is to be expected since this power law is the asymptotic form of equation ( |B8j ) 
for 6 ^> a, a. 
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C. Statistics of W xg 

In this appendix, we describe our treatment of the statistics of W xg taking into account their 
non-gaussian nature and the correlation between measurements at different angular lags. We 
first show how we can approximate the distribution of our data with the multi-variate gaussian 
distribution (§ |C.1| ). We then estimate the covariance of W xg from the data (§ |C.2| ). Finally, in 
\ we apply these results to the fit of our measurement to the Poisson+Clustering model of 



equation (22). While most of the following treatment is standard, we have used the formalism and 



notation given in Lupton (1993). 



C.l. Multi-Variate Gaussian Approximation 

Let us consider the distribution of W xg at a given angular lag 6{. For simplicity, let 
Wfj = W xg j(a p , 0i), where f refers to one of the Nf fields in our data set and i = 1, . . . , N$. Let 
us denote the intrinsic mean and rms standard deviation of the distribution of Wfi as W% and 
(JWi, respectively. 

One complication comes from the fact that Wf , is, in general, not normally distributed 
(see figure |2|). Nevertheless, by the central limit theorem, the distribution of the average 

Nt 

= HidiWf^/Nf is close to a normal distribution with mean W, and standard deviation 
a Wi I \J Nf , if Nf is sufficiently large. In the case of our non-cluster-target data set, Nf is equal to 
33 and should be large enough to satisfy this condition. 

Another complication comes from the fact that the distribution of the Wj's at different i's are, 
in general, not independent. This can be seen in figure ||. This is due to the fact that the same 
fields are used for the measurements at each 9i and that the PSF (and the galaxy/X-r ay-sources 
CCF w gp , if present) extend beyond one pixel. As a result, the joint distribution of the Wi's can 
not be assumed to be simply the product of normal distributions. 

We choose the simplest approximation to the joint distribution of correlated normal variables, 
namely the multi-variate normal distribution. Under this approximation, the joint Probability 
Distribution Function (PDF) for the averages W = (Wi, . . . , Wn b ) t is 

p(W) = . 1 =exp(-(W- W) r V~ 1 (W- W)/2) (CI) 
^{7nr) N 9\V\ 

where W = (WT, . . . , W Ng ) T , and V is the N e x N e covariance matrix which is given by 

V = ((W - W)(W - W) T ). (C2) 

The quantity \V\ is the determinant of V, and the superscript T stands for transpose. An estimate 
for V from the data is derived in the next section. 
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C.2. Estimate for the Covariance 

Let us consider a large number N r of realizations of our experiment. This would involve the 
measurement of W X g(0i) for N r x Nf fields. Our actual measurement is, say, the r = 1 realization. 
Let us denote the resulting measurements as W r ji, and their average within realization r as 
W Tt i = (W r ji)f. In the previous expression and hereafter, a bracket of the form (• • -)i denotes an 
average over the index I. 



We can then express equation (C2) as 



Vij ~ ((Wr,i-W i )(W r , j -W j )) r 

i N f N f 

= ^EEW/ir-Wi. (C3) 
iv / /=l/'=l 

When / ^ /', Wr,/,i an d Wr,/' ,j correspond to measurements in different fields and are thus 
independent. Hence, J2f^f'{W r j,iW r jij) r = Nf(Nf — l)WiWj. On the other hand, when 
/ = /', W r jj and W r fj involve measurements in the same field and are thus not independent. 
To estimate this term from the available data, we approximate (T^r,/,<Wr,/j)r — {Wij^Wijj) f, 
where, as noted above, the r = 1 realization corresponds to our actual experiment. Then, 
Ef=f(W r ,f,iW r , f j) r ~ EfW 1}f ,iW ltf j. 

By combining the f ^ f and f = f terms into equation ( p3[) and after using the 
approximation Wk — (Wf,k)f = Wk, we finally obtain 

Vij ^ ^ ((W U - Wi)(W fJ - Wj)) f , (C4) 

where the r = 1 subscript has been dropped. Note that the diagonal elements, Va are simply equal 
to {(Wfj - Wi) 2 ) f /N f ~ af/Nf, as they should be. 

Our measurement of the covariance matrix V for the non-cluster data set is displayed in 
figure 13. Correlations at low (6 <; 3a p ) and at high (9 ;> 8a p ) angular lags are clearly visible. 



C.3. Fit to Linear Model 

In § |5.3| , we need to fit our measurements of W to the function f(0;a) = 

Qpoisson 

A(0) + 

Ojclustering 

E(0), where a = (a poisson , a c i ustering ) T (Eq. [||]). This corresponds 
to a fit to the N a = 2 parameter linear model (see eg. Lupton 1993, p. 81) 

W = f(0;a) + eEMa + e, (C5) 

where M is a Ng x N a matrix with rows equal to (A(#j), E(0j)), i = 1, . . . ,Ng, and e is the 
A^e-dimensional vector of errors with (e) = and covariance matrix (ee T ) = V. 
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Within the multi-variate approximation of equation (CI), the probability density for our data 
set W to result from a model with parameters a is 

p(W;a)= . 1 e~ x2 ^l\ (C6) 
^(27r) N e\V\ 

where, 

X 2 (W; a) = (W - f (0; a)) T V~ x (W - f (0; a)). (C7) 

The best fit to our model correponds to the value a of a for which p(W; a) is maximum. 
Setting the derivative of X 2 to yields 

a = (M T V _1 M)- 1 M T V'~ 1 W. (C8) 

Note that a is a bilinear combination of multi-variate normal variables and thus itself follows a 
multivariate normal distribution (with two variables). Thus, the PDF of a is 

P(&) = y=j exp (- (a - af IT 1 (a - a)/2) , (C9) 

where U = ((a — a) (a — sL) T ) is the covariance matrix of a, and a corresponds to the true 
distribution of W, i.e. W = f(0;a). By combining equations ( |C8| ) and (|C2|), it is easy to show 
that U = (M^V^M)" 1 . Taken individually, each of the aj's are normally distributed with 
standard deviations simply given by 

<r&i = U u . (CIO) 

These standard deviations are the estimates for the la errors of the best-fit parameters taken 
separately. 

A measure of the goodness of fit is provided by X 2 (W;a). This quantity follows a 
X 2 -distribution with Ng — N a degrees of freedom. Its value must be close to Ng — N a for the fit to 
be acceptable. 

Figure [l(] shows the probability contours for X 2 in the a-plane for our fit to the non-cluster 
data set. The best fit parameters (Eq. ||C8|| ) and their associated uncorrelated la errors (Eq. | |C10f ) 
are also shown. 
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Table 1. Field characteristics and cross-correlation results 



Seq. No. a 


POSS b 


R.A. C 


Dec. c 




b d 


t c 

vexp 


Target f 


{I) s 


(JV) h 


W xg (0) 


W xg (7a p ) 


29 


756 


17h 09m 59s 


( la 


09m 


59s 


33.7 


20053 


4 


2.44 


0.194 


-0.048 


-0.015 


270 


1201 


Olh 16m 29s 


08d 


13m 


59s 


-53.7 


15572 


3 


2.49 


0.072 


0.074 


-0.001 


280 


1563 


12h 24m 59s 


09d 


29m 


59s 


71.2 


41941 


3 


3.67 


0.101 


0.081 


0.019 
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3 
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6490 


2 
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3 


3.25 
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-0.026 
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Table 1 — Continued 



Scq. No. a 


POSS b 


R.A. C 




Dcc. c 




b d 




Target' 


(7)8 


(JV> h 


W wg (fi) 


W xg (7a p ) 


10384 


1421 


14h 


44m 


03s 


07d 


41m 


21s 


56.4 


6548 


3 


3.24 


0.133 


0.011 


-0.015 


10393 


1440 


14h 


26m 


33s 


Old 


30m 


36s 


55.1 


2238 


2 


3.06 


0.110 


0.088 


-0.099 


10437 


1051 


14h 


14m 


13s 


09d 


06m 


38s 


62.8 


9294 


2 


3.78 


0.142 


0.011 


-0.032 


10452 


1283 


02h 


35m 


06s 


Old 


45m 


29s 


-51.2 


16227 


3 


3.28 


0.074 


-0.043 


0.051 


10464 


1119 


15h 


32m 


46s 


23d 


40m 


05s 


53.0 


16905 


2 


3.11 


0.131 


0.008 


-0.017 


10474 


1087 


14h 


44m 


35s 


lid 


47m 


59s 


58.8 


9729 


3 


3.93 


0.145 


0.053 


-0.019 


10533 


1056 


16h 


48m 


41s 


05d 


04m 


59s 


28.9 


33324 


2 


4.39 


0.208 


-0.005 


0.018 


10632 


363 


03h 


23m 


37s 


02d 


14m 


46s 


-42.4 


9269 


1 


4.00 


0.121 


-0.046 


0.033 


10671 


860 


21h 


44m 


01s 


04d 


20m 


30s 


-35.3 


10532 


3 


3.36 


0.114 


-0.066 


0.005 


10682 


1195 


OOh 


14m 


19s 


16d 


19m 


59s 


-45.5 


16733 





2.56 


0.104 


0.018 


-0.030 



^Einstein Observatory sequence number 
b POSS E-plate number 
c Right Ascension and Declination (1950) 
d galactic latitude (°) 

e Einstein-IP 'C Exposure time for VG=l-2 (s) 

f Einstein target type: 1 = galactic sources, 2 = discrete extragalactic sources (galaxies, QSO,etc), 3 = clusters of galaxies 
and superclusters, 4 = deep surveys; one field, labeled 0, has no target type specified in the Einstein database 

g Mean X-ray intensity (10~ 4 cts s^ 1 arcmin~ 2 ) 

Mean galaxy counts (arcmin - 2 ) 
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Table 2. Zero-lag results. The values for W xg (a,0) are given for the real (hard-galaxies) and 
control data sets. The cell size a is shown in units of the pixel size a p = 64", and the errors 

correspond to la. 



a/a v 


hard-gals 


hard-gals(sol) 


soft-gal 


hard-stars 


scrambled 


1 


0.045±0.012 


0.029±0.012 


0.013±0.011 


0.002±0.004 


0.004±0.006 


2 


0.037±0.011 


0.030±0.010 


0.013±0.010 


0.001±0.003 


0.003±0.006 


3 


0.035±0.011 


0.026±0.009 


0.013±0.009 


0.004±0.003 


-0.005±0.005 


4 


0.019±0.008 


0.020±0.009 


0.015±0.006 


0.000±0.002 


0.002±0.004 


5 


0.023±0.009 


0.020±0.007 


0.020±0.007 


0.001±0.002 


0.003±0.004 


6 


0.022±0.008 


0.019±0.006 


0.004±0.005 


0.001±0.002 


0.001±0.003 


7 


0.021±0.007 


0.017±0.006 


0.008±0.004 


-0.001±0.001 


0.002±0.003 


8 


0.020±0.007 


0.016±0.007 


0.010±0.004 


O.OOliO.001 


0.000±0.002 



Table 3. Non-zero lag results. The value for W xg (a p , 6) at different values of 6 is given for the 
whole data set (all) and for the subsamples containing only fields with cluster and non-cluster 
targets. The results for the scrambled set are given in the last column. The angular lag 6 is 

shown in units of the pixel size a p = 64". 



9/a p 


all 


cluster 


non-cluster 


scrambled 





0.045±0.013 


0.064±0.023 


0.029±0.013 


-0.002±0.010 


1 


0.028±0.008 


0.041±0.016 


0.018±0.007 


-0.000±0.005 


2 


0.018±0.009 


0.037±0.016 


0.002±0.007 


-0.000±0.004 


3 


0.013±0.007 


0.022±0.012 


0.005±0.005 


0.001±0.003 


4 


0.021±0.007 


0.042±0.013 


0.002±0.004 


-0.001±0.003 


5 


0.014±0.007 


0.030±0.013 


0.000±0.003 


0.002±0.003 


6 


0.014±0.007 


0.031±0.013 


-0.000±0.005 


0.002±0.003 


7 


0.015±0.006 


0.034±0.011 


-0.001±0.004 


-0.001±0.003 


8 


0.016±0.005 


0.028±0.011 


0.005±0.004 


0.000±0.003 


9 


0.014±0.005 


0.027±0.010 


0.003±0.005 


0.001±0.003 


10 


0.019±0.005 


0.030±0.009 


0.010±0.006 


0.003±0.003 


11 


0.012±0.006 


0.025±0.010 


0.000±0.005 


0.003±0.004 


12 


0.014±0.005 


0.023±0.008 


0.006±0.006 


0.004±0.004 
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Fig. 1. — Example for a XRB galaxy field pair. This field corresponds to Einstein-IP C sequence 
number 2598. The axes give offsets from the center positions with R.A. and Dec. increasing in the 
negative-x and positive-y directions, respectively. The .81-3.5 keV XRB intensity is shown as the 
grey-scale map. Positions of APM galaxies with 13.5 < E < 19.0 are shown as filled circles. 
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Fig. 2. — Distribution of W xg (a p ,0) for the real and scrambled data sets. The mean values are 
.04 ± .01 and —.001 ± .008 for the real and scrambled sets, respectively. 
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Fig. 3. — Zero-lag results. The values for W xg (a, 0) is shown as a function of cell size a for the real 
(hard-galaxies) data set and for the control data sets (hard XRB with high solar contamination and 
galaxies, soft XRB and galaxies, hard XRB and stars, hard XRB with galaxies for scrambled field 
pairs for which the dashed lines delimit the ±la uncertainties). The dotted-dashed line correspond 
to the zero-lag values expected from our phenomenological fit to the non-zero-lag results (see text). 
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Fig. 4. — W xg (a p ,0) for the real data set vs. various quantities: (a) Einstein-IP C exposure time; 
(b) absolute Galactic latitude; (c) mean XRB flux in the .81-3.5 keV band (VG=l-2); and (d) mean 
number of APM galaxies with 13.5 < E < 19. 
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Fig. 4b.— 
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Fig. 4c— 
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Fig. 4d.— 
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Fig. 5. — Non-zero lag results. W xg (a p ,6) vs. for the real and the scrambled data set. The 
dot-dashed curve shows the result of a phenomenological fit used for a comparison with the zero- 
lag results. The dotted curve corresponds to the best fit to the Poisson + clustering terms, while 
constraining the fit parameters to be positive. 
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Fig. 6. — Illustration of the correlation between measurements of W xg (a p ,9) at different values of 
9. For each of the 62 fields, W xg at 9 = 10a p is plotted against that at 9 = lla p . 
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Fig. 7. — Field with sequence number 8926. For clarity, the .81-3.5 keV XRB intensity (a) and the 
distribution of galaxies with 13.5 < E < 19.0 (b) are shown separately. This field has the highest 
non-zero lag correlation in our data set. A cluster of galaxies is visible in the upper left-hand corner. 
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Fig. 7b.- 
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Fig. 8. — Zero-lag vs. Non-zero lag values of W xg for each field. W xg (a p , 7a p ) is taken as 
representative of the non-zero lag value. Fields with clusters and non-cluster Einstein targets 
are displayed as stars and diamonds, respectively. 
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Fig. 9. — Non-zero lag results for cluster and non-cluster targets separately. The value of W xg for 
all fields combined is redisplayed for comparison. The solid line is the best fit of the Poisson + 
clustering terms to the non-cluster target field values. 
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Fig. 10. — Confidence contours for the fit of the Poisson + clustering terms to the non-cluster fields 
values of W xg (a p ,6). The best fit parameters and their associated error bars are displayed as the 
cross. 



~ i 1 1 r 



~ i 1 r 



10" 



E 
o 

I 

w 10" 

cn 

i_ 

>' 
CD 

O 

^ 10" 
l 

CN 



10" 



□ E + SO 
10 )K E + SO (upper limits) 
+ W xg correlation 



□ 



□ 



F 



□ □ 











□ 




□ 























□ 




_l I I L 



_1 I L_ 



20 



15 



10 



B 



Fig. 11. — Apparent B-magnitudes and X-ray fluxes for the galaxies in the catalog of Fabbiano et al. 
(1992). The values are shown separately for: (a) Spirals and Irregulars, and (b) Ellipticals and SOs. 
The dashed lines correspond to the la limits for each of these types of galaxies, assuming a constant 
X-ray-to-optical flux ratio. The value of the mean galaxy flux derived from our cross-correlation 
analysis is indicated as the cross. 
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Fig. 12. — Functions used in the interpretation of W xg for parameters relevant to our study, i.e. 
for a = a p = 64", a = a ps f = .5', and (3 = .7. The dot-dashed curve is the PSF. The solid curve 
is A(a,a,8). The crosses are the result of the Monte-Carlo evaluation of S(a, a, (3, 0) with error 
bars corresponding to 2a limits. The dashed curve is the analytic function used to model H. For 
comparison, the asymptotic form (#/£o) -/3 is displayed as the dotted line. All functions have been 
normalized to one at 9 = 0. The actual normalizations are A(0) ~ .551 and 2(0) ~ 1.480. 
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